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I.  INTRODUCTION 


The  heat  transfer  occurring  at  the  surface  of  a  solid  in  contact 
with  a  fluid  at  different  temperature  is  frequently  modeled  by  the 
radiation.  Newton  cooling  or  convective  heat  transfer  condition 

3u 

K  iTP  '  H  (Uf  •  us);  (1>1) 

here  IT  is  the  inward  normal  from  the  surface  of  the  solid,  Uf  and  u5  are 
the  temperatures  of  the  fluid  and  solid  respectively  at  their  interface, 
k  is  the  coefficient  of  thermal  conductivity  for  the  metal  and  H  is 
the  effective  heat  transfer  coefficient  (which  might  include  radiation 
terms).  In  most  cases  H,  K  and  Uf  are  not  known  a  prior:  since  they 
may  vary  with  existing  conditions  and  their  determination  becomes  a 
major  portion  of  the  problem.  However,  to  make  progress  and  to  under¬ 
stand  the  role  of  these  parameters  more  completely  we  may  idealize  the 
problem  by  assuming  that  h  =  H/K  is  a  known  positive  constant  and  up  is 
a  given  function  of  time.  Ultimately  these  are  not  serious  restrictions 
when  incorporated  into  a  larger  scale  fluid  flow  calculation  since  h  and 
Uf  can  be  reset  over  each  sufficiently  small  time  increment  of  the  corn- 
put  at  ion . 

The  objective  of  the  present  report  is  the  construction  of  a  number 
of  exact  solutions  for  the  heat  or  diffusion  equation  which  satisfy 
boundary  condition  (1.1).  These  functions  will  be  used  subsequently  in 
a  separate  report  to  develop  an  effective  numerical  algorithm  which 
models  the  heat  transfer  occurring  at  the  bore  surface  of  gun  barrels. 

In  this  report  however  we  shall  not  be  concerned  with  engineering  appli¬ 
cations  but  formulate  the  problem  as  a  purely  mathematical  one.  Our 
analysis  will  suppose  that  the  temperature  of  the  solid  depends  on  a 
single  spatial  coordinate  x  and  on  time  t  and  that  the  solid  exists  in 
the  region  x  >  0.  Denoting  the  temperature  of  the  solid  by  u(x,t)  the 
boundary  condition  (1.1)  may  be  rewritten  in  the  form 


I-h  u  (O.t)  =  uf  (t) 


where  1. 

n 


is  the  differential  operator  defined  by 


(1.2) 


I-h  U  *  U 


1  du 
h  dx 


(1.3) 


with  h  «  H/K  >  0. 
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In  Section  11  we  shall  define  our  problem  precisely  and  introduce 
the  transformation  which  is  the  inverse  of  1.^.  In  Section  111  we 
find  transformations  for  the  functions  and  Mv*  which  were  introduced 
in  previous  work1**.  Finally  in  Section  IV  explicit  solutions  of  the 
problem  are  obtained  for  specific  choices  of  the  boundary  and  initial 
data.  We  shall  make  frequent  use  of  the  results  in  Reference  2  and  shall 
assume  that  this  is  available  to  the  reader 


II.  ('ROB  1.1. >-1  STATEMF  VI  AND  1111  TRANSFORMATION  1. 

h 

Consider  the  following  problem  defined  in  the  quarter  plane  x  >  0, 

t  >  0. 


du 

a  • a 

3  u 

3x“ 

x  >  0,  t  »  0 

(2.1) 

u  (x ,0)  -  f 

(*) 

x  >  0 

(2.2) 

1^  u  (0 ,  t )  *  t. 

(t) 

t  »  0 

(2.3) 

where  a 

•  0  and  I,  is  defined 

by  (1.3). 

Ti 


To  obtain  solutions  for  this  problem  we  shall  proceed  using  purely 
formal  derivations.  That  is,  we  shall  not  be  concerned  with  mathematical 
legitimacy  hut  presume  that  all  functions  in  question  exist  and  are 
sufficiently  smooth  to  justify  whatever  operations  are  performed. 
Vindication  of  this  procedure  follows  a  posterior:  by  noting  that  the 
functions  thereby  obtained  do  in  fact  satisfy  the  desired  equations. 
i>n  the  other  hand,  the  important  questions  of  uniqueness  and  continuous 
dependence  on  the  data  will  not  be  considered,  in  fact,  the  solution  to 
problem  i2  1)  -  (2.3)  is  not  unique  unless  an  additional  exponential 
growth  restriction  is  stipulated.  (Consider  in  this  regard  the  example 
given  by  Friedman*).  Such  restrictions  serve  mainly  to  rule  out 
physically  meaningless  solutions  and  will  certainly  be  satisfied  by  the 
functions  constructed  below  which  are  either  hounded  or  grow  more  slowly 
than  a  polynomial  in  x  (as  x  *  »  with  t  held  finite). 

F.  Folk,  "Asy^totic  Expansions  for  the  Solutions  f  Parabolic 
,  fferential  Imitations  oi  th  a  ->w  1 1  Ptirar*- 1  e  r ,  "  rh.l  iiss>  rtati  >»;, 
rtru-nt  Hxthe'iiti  v»,  'huiversi  ty  of  N  '.t.'w  ,  ,V.-:,wrJ ,  ;»  ht .'ar>  , 

Jim 

■.'ohn  F.  Folk,  "Special  function  Solutions  of  The  Diffusion  Equation," 
Technical  Heport ,  B>xl  listic  Eesear\-h  laboratory, 

'Abner  Frio  inrr:,  "Partial  Differential  Equation  of  Parabolic  Type," 
Prentice-Hall  lorp,  Fn^lstoooi  I'ress,  N.J.,  p.  31,  19F4. 


To  solve  problem  (2.1)  -  (2.3)  note  that  the  related  function 


satisfies  the  following  equations 


3v 

"it 


a  v 

a  *T 
3*‘ 


v ( x . 0 )  *  1^  f  (x) 
v(0,t)  *  g  (t) 


V 


x  »  0 ,  t  >  0 

x  >  0 
t  >  0  . 


(2.4) 

(2.5) 

(2.b) 

(2.7) 


Using  linearity  we  may  express  the  solution  to  this  problem  in  the  form 


where  v}-  and  v  are  the  solutions  to  (2.5)  -  (2.7)  with  g  (t)  *  0  and 
f  x  *  0  respectively.  But  these  problems  have  the  well  known 
solut ions 


v  ,(x.t)  ■  /  (H  j  (x-s.  at)  -  H  j  (x*s.  at)][l.^  f)  (s)  ds  (2.8) 

o 

and 

v  (x.t)  •  -Ta  ''  H  _  (x.a  (t-r))  g(r)  dr  (2.9) 

K  o  '* 

where 

H  j  (x.t)  •  (4nt)  1  ‘  exp  (-x‘/4t)  (2.10) 

H  ,  (x.t)  •  .t-  H  (x.t)  *  (-x/2t)  H  ,  (x.t)  .  (2.11) 

The  functions  H  (x.t).  for  all  >  €  R,  were  investigated  in  References 

1  and  2  and  will  be  briefly  reviewed  in  the  next  section.  The  expres¬ 
sion  for  Vf  is  somewhat  unsatisfactory  in  that  it  requires  f(x)  to  be 
differentiable,  although  this  was  not  needed  in  the  formulation  of 
problem  (2.1)  -  (2.3).  *e  can  eliminate  this  difficulty  using  integra¬ 
tion  by  parts,  as  follows 


V-  (x,t)  *  /‘[H  (x-s,  at)  -  H  (x*s.  at)]  f(s)  ds 

t  o  * 1  ' 1 

}  r  (H  ,  (x-s,  at)  -  H  .  l>*s.  at)]  f'(s)  ds 
h  '  - 1  "i 

o 

=  /*  (H  {  (x-;,  at)  -  H  1  (x*s,  at)]  f(s)  ds 

o 

-  1  “  (H  ,  (x-s,  at)  ♦  H  (x*s,  at)]  f(s)  ds 

n  -  * 

o 

no 

♦  I H  j  (x-s,  at)  -  H  j  (x*s,  at)]  f(s)  ls  ,  Q  • 


The  last  term  vanishes  because  of  the  exponential  deca>  of  H  j(x,t) 
prov i Jed 

f ( x I  <  const,  exp  [a  x‘] 

f.  ;•  some  a  •  1  / 4 a t  .  The  remaining  terms  can  be  recombined  in  the  form 

v  ix.t)  =  [1^  H  j]  (x-s,  at)  -  [l-h  H  j  ]  (»*S  at) 

°  (2.12) 


I  x*s ,  at  I )  f (s)  Js . 


The  solution  to  problem  (2.11  -  (2.3)  can  now  be  obtained  by  inverting 

the  transformation  u  •  v  -  L  u.  This  amounts  to  solving  an  ordinary 

n 

differential  equation  and  can  be  accomplished  using  the  integrating 
factor  method ■  The  appropriate  transformation  is  denoted  by  T^  and  de¬ 
fined  for  a  function  w  =  w(x)  by 

IT.  w]  (x)  h  exp  |h  (x-s)]  w(s)  ds.  (2.13) 

h  x 

This  is  well  defined  and  finite  for  all  w(x)  in  the  class  of  functions 
whivh  are  defined  and  continuous  for  x  >  o  and  bounded  by 

w(x)i  <  const,  exp  (h  x) 

for  some  h  h  and  sufficiently  large  x.  The  transformed  function  is 
differentiable  and  can  be  similarly  bounded.  Thus  T^  w  €  whenever 

w  (  s  .  Moreover  we  have 

n 
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Th  w(x) ;  <  h  M(x)  /  exp  [h  (x-s) |  ds 
x 

<  M(x) 

where 

M(x)  *  sup  I  w(s) | :  s  >  x) . 

It  follows  that  if  w  is  hounded  or  vanishes  asymptotically  as  some 
negative  power  of  x  as  x  •  »  then  w  is  also  bounded  or  vanishes  at 
least  as  rapidly. 

Applying  Leibnitz's  rule  to  (2.13)  we  have 

[T^  w ) * ( x )  «  h*  /  exp  (h(x-s) ]  w(s)  ds  -  h  w(x) 
x 

«  h  (Th  w  .  w]  (X) 

and  consequent  1 > 

Lh|Thw].w  (2.14) 

whenever  w  is  ot  class  S.  .  The  transformations  L  and  T.  are  thus  seen 

h  n  h 

to  be  inverses.  If  w(x)  has  a  derivative  w‘(x)  which  is  also  of  class 
"Y  then  using  intergration  by  parts  we  have 

[Th  “I  ( x )  *  h  exp  [hlx-sl]  w*(s)  ds 
x 

■  h  exp  (h(x-s)  ]  w(s) 

s*x 

♦  h*  exp  [h(x-s)]  wls)  ds 
x 

»  h  [Th  w  -  w]  (x) . 

Thus 

|Th  w')  (x)  -  |Th  w)'(x).  (2.15) 

That  is.  the  derivative  of  the  transform  is  the  transform  of  the 
derivative,  this  permits  us  to  reverse  the  order  of  the  operations  in 
2.14>  for  functions  with  derivatives  of  class  S^: 
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Th  lLh  -1  *  Th  -  -  h  Th  l-'l 


Th  w  -  i  iTh  -l 


hi  1Th 


or 


Th  I'-h  “i 


(2.14) 


In  general  if  the  n-th  derivative  w'n  (x)  exists  and  is  of  class  then 
by  repeated  applications  of  (2.14)  we  obtain 


IT,  “ln  ]  (x)  =  (Th  w]ln'  (x) 


(2.1b) 


for  n  ■  0,1,2 . 

he  can  now  return  to  the  original  problem.  The  functions  v,  and  Vi 
are  twice  differentiable  with  respect  to  x,  for  t  •  0,  since  they  are 
solutions  of  the  heat  equation.  Furthermore  they  can  easily  be  shown 
to  be  of  class  S,  (for  each  fixed  t  >  0)  provided  f  is  of  class  and 
g  is  bounded.  Tnus  their  transforms  are  well  defined  and  the  solution 
to  problem  (2. It  -  (2.3)  can  be  written  as 

u  *  ur  ♦  u  ■  T.  v,  4  T.  v  (2.1?) 

f  g  h  f  h  r 

where  u(  and  u  are  the  solutions  of  the  problem  when  R(t)  -  0  and 
fix'  -  0  respectively.  But  applying  to  the  expressions  (2.9)  and 

2.121  and  using  Fubini's  Theorem  to  interchange  the  order  of  the 
integration  and  transformation  operators  yields 


u^(x.t)  -  j  [11  ,(x-s,a  t)  -  H  ,(x*s,a  t) 
o' 


h  Th 


( x ♦ s , a  t))  f(s)  ds 

(2.18) 


u8(,.t) 


Th " 


(x.a  (t-r))  g(r)  dr  . 


(2.19) 


Transformation  of  the  functions  H  ,  and  H  ,  is  permissible  since 

thev  are  not  only  continuous  but  entire  with  respect  to  x  for  t  >  0 
and  vanish  exponentially  as  x'  -*  » .  We  also  note  that  the  identity 
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Th  [w  (**y) )  ■  (Th  w] (x*y) 


(2.20) 


tar  any  function  w  of  class  follows  directly  from  the  definition  of 

T  .  Thus  the  application  of  T.  to  M  ,  in  (2.18)  is  unambiguous  in  that 
h  n  -i 

the  shifting  (by  si  can  be  performed  before  or  after  the  transformation 
(with  respect  to  x) . 

In  the  next  two  sections  we  shall  obtain  explicit  expressions  for 
the  transformations  of  the  functions  14^  (x,y),  i  €  R  and  solutions  to 
problem  12.1)  -  (2.3)  for  specific  choices  of  f(x)  and  g(t). 

III.  T.  TRANSFORMAT I ONS  01  H  ,  H  *  and  v 
h  y  >  n 

let  us  briefly  review  the  salient  properties  of  the  functions  H  , 

(I  *  and  which  were  discussed  in  References  1  and  2.  These  were 
Y  n 

defined  by 

H  (x.t)  =  i  4-w  t ) 1  “  /  v  exp  (-  lx-s)‘/4t)  ds  (3.1) 

Y  ‘  >  ! 

O 

for  >  •  -  1  anJ  recursively 

H  (x.t)  •  ~  H  ,  (x.t)  (3.2) 

1  >X 

for  i  s»  -  1.  (The  latter  equation  also  holds  for  >  ■»  -  1 )  .  In 
part icular 

H  (x.t)  »  (1/2)  erfc  (-  x/2  /T) 

o 

H  j  ix. t)  *  (4-t)  1  ‘  exp  (•  x*/4t )  (3.3) 

H  ,  (x.t)  «  (-  x/2t)  H  j  (x,t) 

where 

erfc  (i)  =  (2//«  )  /  exp  (-  s~)  ds  . 

y 

!he  associated  functions  H  *  were  given  by 

H  *  (x,t)  ■  (-  x.t)  (3.4) 
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and  for  any  integer  n  the  heat  polynomials  were  defined  by 

vn  (x.t)  =  (x.t)  ♦  (-1)"  »n  (x.t);  (3.5) 

for  negative  integers  these  vanished  in  view  of  the  identity 

H*  (x.O  *  (-1)"*1  Hn  (x.t)  (3.6) 

tor  n  *  -l,  -2 .  All  of  these  functions  were  seen  to  be  infinitely 

differentiable  in  x  .ind  t  for  t  ■*  0,  and  satisfy  the  formulas 


f-  H  *11  , 

3x  >  Y-l 


H  -  H  , 
3t  >  1-J 


•  • 

H  *  -  H 
3x  >  >-l 


3  •  • 

—  n  =n  , 
at  y  y -2 


-  v  *  V  . 

ax  n  n  - 1 


- —  v  *  v 
3 1  n  n  -  2 


Thus  the\  are  all  solutions  of  the  heat  or  diffusion  equation 


3  u  3u 
.  :  ’  it 


(3. 10) 


t  r  t  •  and  satisfy  the  initial  and  boundary  values 


H  (x,0)  »  h  (x 

Y  Y 


H.  ( 0 , t )  *  *  h  /:(tl  (3.11) 


M  ( x  ,0)  *  h  ( x ] 
Y  Y 


v  ( x  , 0 )  =  x  /n! 
n 


(O.t)  *  ^  h.  f2  (t)  (3.12) 


,  (t )  n  even 


"  I  0 


n  odd 


(3.13) 


where  h  and  h  denote  the  jump  functions 


\<») 


I 


x  »  0 
x  v  0 


1 


h  ( x )  ■  h  t  -  *  ) 

Y  Y 


(The  reciprocal  of  the  factorial  function  is  defined  tor  all  ,6  K  by 


>!  r ( > • 1 ) 

where  denotes  the  usual  f.awna  function.) 

from  the  discussion  in  Reference  1  it  is  known  that  for  each  i  C  R 
there  exist  a  constant  R  ■'  such  that 


H  x,t  *»  k  [h  ( x)  ♦  /t  !  . 

Y  y  Y 


(5.14 


rhus  it  is  clear  that  the  functions  H  x.t  ,  H  'X.t  and  v  x.t  are 

Y  Y  n 

-t  -  lass  with  respect  t  o  x  ,  for  each  fixed  f  >  0,  and  their  trans- 

:  -r~.it  .  :is  are  well  defined  for  convenience  we  shall  use  the  notation 


T.  " 

h  y 


(3.15 


2  •  T.  H  (3.  lh 

Y  h  Y 

:  e  here  that  in  contrast  ti  the  simple  formulas  2.1<>!  and  (2.20)  tl 
definition  of  T.  2.13)  implies  sn  general  that 

Th  ‘-*'1  *  lTh  ('X)' 

Thus  the  t  ransf.-rmat  ion  of  M  is  not  Riven  hv  2  (x.t  I  -  2  (-x.t).) 

Y  -  Y  Y 

Using  differentiation  formulas  (3.21  -  (3.9)  and  the  hound  (3.14  to 
w  uniform  convergence  of  the  defining  integral  we  arrive  at  the 
f  -  - rmu 1  a s 


Jx  “y  “y-1 


at  ">  y-2 


for  t  These,  in  turn,  imply  that  2  and  2  are  solutions  of  the 


heat  equation  1.3.10)  for  any  >  e  K.  Applying  T^  to  the  initial  values 
t  3.  11)  and  (3.  12)  we  have 


2  ( * , 0)  ■  h  C  eh(X,S)  *1  ds  (3.19) 

>  *  >  < 

X 

and 

:*  ( x  .0 )  -  0  (3.20) 

tor  \  •  0.  The  former  expression  can  he  explicitly  evaluated  when  i  is 
.»n  integer  by 


n  -  0,1.2 


n  •  -1.-2,  .  .. 


(3.21) 


It  would  he  useful  to  obtain  explicit  functional  form-,  for  these 
tr.insforn.it  inns  and  for  integer  values  of  >  this  is  possible,  .is  we  now 
show  first  consider  i  ■  -1.  which  may  he  regarded  as  a  fundamental 
case,  since  H  ,  is  the  fundamental  solution  of  the  heat  equation.  The 

functions  )l  ,  and  H* (  are  identical  and  their  transformation  is  given  hs 


j  (x.t)  »  1  j  (x.t) 


er  fs 


m 


exp  (hx  •  h  t )  (322 


or,  alt ernat i v e 1 s , 


-  1 


(x.t) 


h  H 


- 1 


(x.t) 


M*  ( x ♦ 2ht  .  t) 

o  _ 

M  'j  (  x  *2Vit  ,  f) 


(3.231 


for  t  >  0.  To  verifs  this  note  that  the  former  expression  mas  he 
differentiated  to  give 


;  ( x  .  t )  *  t-h  .4-~t)  exp  [  -  (,x*-ht  )  *■  /  4 1  ]  exp  |hx»h‘t 

*  x  -  1 


♦  \  erfc  |  (x».’ht  )/«'4t  ]  exp  [hx*h*t] 


( - h / • 4 * t )  exp  [-(xW/4t)]  ♦  h  Z  .  (x.t) 


-h  [H  t  -  Z  j ]  (x.t) 


Thus  the  function  Z  ^  defined  by  (3.ZZ)  clearly  satisfies 


S  =-l  ’  H-1 


V  te  furthermore  that  VZZ;  may  he  written  as 


x.t)  •  ;  erfc  (:)  exp  exp  ( -  x‘/»t ) /.  4-t 


where 


(  x  • h t  »4t  . 


lor  1  a rue  i  the  standard  asymptotic  formula 


then  implies  that  as  x  *  ■*  with  t  '  0  fixed  we  have 


i  '  '  -  <  x  ^3 1 ) 


which  vanishes  exponentially.  It  follows  that  the  expression  defined 
hv  t.ZZ  is  of  class  and  must  indeed  he  the  transform  of  H  j. 

_  • 

»e  are  now  able  to  determine  functional  forms  for  Z_  and  Z  for 

n  n 

’thrr  lntc^rr  valuer  of  n  with  the*  help  of  a  useful  recurrence  relation. 
Using  (3  '  and  the  definition  of  1^  we  have 


(3.24) 

This  permits  us  to  express  in  terms  of  I  j  and  the  functions  : 

»  between  -  1  and  n : 


-3  . 


1 


h“'  H  ♦  h'1  H.  •  H, 


•t  ’  »•“  2.,  *  h'  Ho  *  Hl 


:  -  h'1  :  ,♦  h 

o  - 1  o 


■-1  -1 


(3.25) 


h  2.,  *  h  H-1 


j  '  »>*  2.,  *  h‘  f'-i  -  h  H  - 2 


2.4  ’  h'  :-l  *  h'  H-1  *  h'  H- 2  *  h'  H-3' 
In  general,  using  induction,  we  have 


-n- 1 


n 


’-1 


k*o 

0 

-n- 1 


hV'n  Kv 


\  v-i 


h"n'V  H 


-k 


n  >  0 

n  ■  -1 
n  <  -2 


(3.26) 
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Thus  , 


T.  v  -  T.  |ll  ♦  C- 1  )n  H* ] 
h  n  h  1  n  nJ 

.  n  L 

. -n-l  ,  r  ,k-n  ,, 

.  h  :  •  ;  h  Hk 

k«0 

-  h'n_1  :  t  ♦  i-nK  hkl'n  ii* 

•  l  |M.  ♦  t- 1  )K  »£]. 

k*o 

n  , 

T.v  -  l  hK  n  v.  .  (3.30) 

h  n  V 

k«o 


\  Iced  t  >rm  <•»,  r<  ■  •  i  n  such  as  ,  3.22'  has  not  hern  obtained  for 
r  1 1  »hen  i  i  not  .in  integer.  However  an  alternative  series 

f 

repre sent  at i on  for  Z  can  be  used  in  such  cases: 


* 

v 


(-h)K  H*#k  (k.t). 


1  is  i  ng  i  nequa  1 1 T  s  3.14  it  is  easy  to  show  that  for  any  v  €  R  this  series 
is  umf  ralv  convergent  in  compact  subsets  of  *  •  0,  t  >  0.  Thus  from 
3 .  ?  >.e  see  that  term  bv  terra  differentiation  is  possible  and  that  the 

eric  satisfies  formulas  3  1 " i ,  the  heat  equation  (3.10)  and 


■  H* 

T 


for  each  i  c  R.  Along  x  -  1  this  series  takes  on  the  values 


whereas  along 


:  (O.t) 

y 


m 

I 

k-1 


((>•1072)! 


0  it  vanishes  for  *  >  0.  as  it  must  according  to  (3.20). 


IS 


Fhc  analogous  series  for 


would  seem  to  be 


Y 

l  «  k  U.t) 

k»  1 

which  converges  and  satisfies  the  appropriate  equations.  However,  using  (3.11) 
this  series  can  be  shown  to  grow  exponentially  as  x  «->  and  thus  cannot  be  of 
class  5^  ,  even  though  the  individual  terms  belong. 

in  conclusion  we  note  that  although  our  interest  has  been  in  the 
region  x  >  0,  t  >  0  the  functions  Z  and  Z  can  be  defined  in  the  domain 

x  '  0,  t  *  0  as  well,  and  in  some  cases  along  t  ■  0  also. 

IV.  SOLUTION  01  BOUNDARY  VALUE-  PROBLEM 

let  us  now  apply  the  results  of  Section  111  to  problem  (Z.l)  -  (2.3). 
Substituting  VIS'  into  (Z.l'i  and  (2.18)we  have 


u  *  u 

f  ‘  “g 

(4.1) 

with 

u,  x ,  t  (H  Ix-s.at 

-  » 

o 

y 

-  H  (x*S ,at)  -  Z 

,(»*'.«t  1 

f  (s 

)  ds  (4.2) 

U  (  *  .  t  i 

ft 

,  rt  - 
•  -  2a  * 

2 ( x . a ( t -r) 

)  g( r)  dr 

(4.3) 

0 

where  Z  is  given  by  ( 

•  i 

22 )  and 

•  •>  1 

These  formulas  represent 

the  general  solution 

for  problem 

(2. 

1)  -  (2.3). 

Next  we  ma\  consider 

T  he  *|-*ei  l  f 

u  choices 

of  boundary 

and 

initial 

data 

f  (  X  )  • 

h  ( x- x  ) 

Y  o 

(4.4) 

R(t)  • 

\/2(t) 

(4.5) 
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Hu-  solutions  ur  and  u  for  these  cases  can  be  obtained  more*  easilv  bv 
f  K 

first  solving  for  v^.  •  1^  u^.  ang  vt  ■  u  and  transforming  by  than 

by  using  formulas  (4.’)  and  (.4 .  .^ )  directly.  These  functions  are  required 
to  satisfy  Equation  (2.S)  with 


v-(x.O)  ■  h  ( x * x  )  -  j-  h  .  ( x - x  ) 
t  Y  o  h  >  - 1  o 


vf (0 ,t ) 


and 


V  (x  .01 

K 


V°'t) 


•  t 


7,7m 


hut  in  previous  work.-  »e  have  seen  that  these  problems  have  the  solutions 


V.(X.t)  -  M  (x-x  .  at)  -  H  ( x • x  ,  at 
f  Y  O  Y  O 


1  1  * 

.  M  ,(x-x  ,  at)  ♦  — H  ,(x*x  .  at) 

h  Y-J  o  h  y - 1  o 

and 

v  (x.t)  .  2H* (x ,  at ) . 

K  Y 

These  can  he  directly  verified  using  iVll  and  (VI.1)  Bv  applying 

the  transt  rm  T.  to  v  .  and  y  we  obtain 
h  f  g 


'  i:.  h  at  i:>  h i; 


•  Il*h  }  -  "*h 


i  2  .(X'x  .it) 

n  y  - 1  o 


or 


u  (x.t)  ■  H  (x-x  ,at)  -  H  ( x »x  .at)  •  c  I  ,(x«x  ,at) 
t  Y  O  YO  h  Y*  *  O 


M.b) 
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u  (*,  t  ) 

g 


(4.7) 


and 


2  :t(x.at). 


*e  iiin  illustrate  our  procedures  by  giving  as  examples  the  solutions 
t  the  following  two  important  particular  choices  for  f(x)  and  g(t): 

f  ( x )  ■  1  f  ( x )  ■  0 

g(t)  •  o  and  g(t)  -  1. 

1 r  "i  l  t>  and  1  '  these  are  seen  to  have  the  respective  solutions 

U-lx.t)  ■  H  x .  a  t )  •  H  *(x,at)  ♦  r  Z  ,(x,at) 
i  o  o  n  - 1 

* 

u  x.ti  *  2  Z  ( x  ,  a  t ) . 
g  o 

sing  \.22  and  t .  to  expand  these  yields 

1  -  1  , _  ? 

i.  *  .  t  .  erfc  »  » 4a t  -  -  erf  c ( x  >4at)  ♦  r-  I  , (x .at ) 

.  .  h  - 1 


e  r  t  (  x  •  4  a  t 


u  «. .  t  e  r  f  c  :  x  / »  4  a  0  -  erfc 

g 


'  ' 

(» ♦  .’hat  \ 
/4a?  / 


xp ( hx  ♦  h* at 


explhx  ♦  h'at). 


1  •  t  *.  i  ■  !.-m>  these  lut  i  ■  ns  ire  seen  to  he  identical  with  the  results 
of  Cars  law  and  1aeger“. 


\  CONCLUSION 

The  transformation  defined  hi  .'.It1  was  introduced  in  order  to 
■overt  '■  lut  i  ns  >f  the  heat  equation  with  known  boundary  values  into 
-'lilt  ions  itisfiing  the  onvective  heat  transfer  condition  (1.2)  along 

»  In  particular  the  functions  M  and  H  were  transformed  into 

y  y 

r  ;;  ti  •  and  Z  '  general  s.  lut  ion  4  1  -  4  t  was  found  for  the 

Y  Y 

mve,  t  ive  heat  transfer  problem  (2.1  -  12.'  and  specific  solutions 

were  found  for  particular  choices  4.4  and  i4.Sl  of  the  boundary  and 
i m  1 1 a  1  data . 


1  .  .  iraZ  nc  mi  .  .  \  J wjer,  ”  '  mht.'ti  -n  of  He.it  in  Solute,”  Oxford 

'Usmndon  Pme a,  19i>9,  PnJ  Addition  (pp  ?0-?2). 
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!ne  basic  problem  l».l)  -  ( J . 3 )  studied  in  this  report  was  posed 
tor  the  semi-infinite  domain  *  >  0,  t  >  0.  However  all  of  the  functions 
encountered  were  extendable  into  x  <  0,  t  >  0  and  thus  can  be  used  in  a 
more  general  context.  They  can  also  be  used  to  develop  approximate 
solutions  for  more  general  equations  in  bounded  spatial  domains.  In  a 
subsequent  report  we  shall  use  these  functions  to  construct  a  numerical 
algorithm  which  models  the  heat  transfer  occuring  at  the  bore  surface 
of  a  gun  barrel. 
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